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Abstract 

The present paper is an extension of a recent work (Bhattacharya et al. 2010) 
to the Einstein-Strauss vacuole model with a cosmological constant, where we 
work out the light deflection by considering perturbations up to order M 3 and 
confirm the light bending obtained previously in their vacuole model by Ishak ct 
al. (2008). We also obtain another local coupling term __I_M___ related to A, in 
addition to the one obtained by Sereno (2008, 2009). We argue that the vacuole 
method for light deflection is exclusively suited to cases where the cosmological 
constant A disappears from the path equation. However, the original Rindler- 
Ishak method (2007) still applies even if a certain parameter 7 of Weyl gravity 
does not disappear. Here, using an alternative prescription, we obtain the known 
term — ^M, as well as another new local term 37r I M between M and 7. Physical 
implications are compared, where we argue that the repulsive term — can be 
masked by the Schwarzschild term 2M. m the halo regime supporting attractive 
property of the dark matter. 



I. Introduction 

Recently, we confirmed the Rindler-Ishak method (Rindler & Ishak, 2007) 
by calculating light bending in a more general solution, viz., the Mannheim- 
Kazanas-de Sitter solution of Weyl conformal gravity (Bhattacharya et al. 2010) 
that contains two parameters A and 7, the latter is assumed to play a promi- 
nent role in the galactic halo populated by dark matter. The method indeed 
delivered the effect of 7 exactly as it has been known in the literature for long. 
Subsequently, the A— effect has been calculated by Ishak et al. (2008) within 
the framework of Einstein-Strauss vacuole (its earlier incarnation is the Kottlcr 
vacuole). Our broad aim here is to examine how the effects of both A and 7 
appear from suitable considerations in the Weyl conformal gravity. 

Galaxies or clusters of galaxies (hereinafter called lenses for brevity) have 
been modelled as residing in the womb of a de Sitter vacuole much larger than 
their sizes (Ishak et al. 2008). The model assumes that, for a given lens, the 
boundary radius r b of the vacuole is determined where the spacetime transitions 
from a Schwarzschild-de Sitter (SdS) spacetime to a cosmological Friedman- 
Robertson- Walker (FRW) background. Further, all the light-bending occurs in 
the SdS vacuole and that once the light transitions out of the vacuole and into 
FRW spacetime, all bending stops. Ishak et al. (2008) showed that the effect 
of the cosmological constant A appears inside the vacuole in the bending of 
light by different lens systems. They also obtained an upper bound on A, using 
observational uncertainties in the measurement of the bending of light, which 
turned out to be only two orders of magnitude away from the cosmologically 
determined value. For a lens of mass M and radius R, Ishak et al. (2008) 
obtained light deflection up to second order in together with a A— repulsion 
term (= -ARr b /6). 

The purpose of the present article is to confirm the light deflection in the 
second order by using perturbations up to third ordeiQ- When we do that 
we come up with new extra terms, while confirming the most interesting term 
— Ai?rfc/6. Additionally, one might like to have an idea of how the presence of a 
conformal parameter 7 would affect light deflection by the lenses. To this end, we 
use an alternatuive orescription for the more general exact Mannhcim-Kazanas- 
dc Sitter solution of Weyl conformal gravity that includes the parameter 7. Pure 
SdS vacuole with only M and A is readily recovered at 7 = 0. New coupling 
terms arising out of the invariant angle have been obtained. 

II. The solution and the approximation scheme 

One well discussed solution that contains both the conformal 7 and dS 
A effects is the Mannheim-Kazanas-de Sitter (MKdS) solution (Mannheim & 
Kazanas 1989; Mannheim 1997, 2006) of Weyl conformal gravity field equa- 
tions. The metric is given, in units G = cq = 1, by (see e.g., Edery & Paranjape 

1 We are indebted to an anonymous referee for suggesting that the correct deflection follows 
from (at least) third order calculation. 
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1998): 

dr 2 = B(r)dt 2 - A{r)dr 2 - r 2 {d6 2 + sin 2 9dcp 2 ), (1) 

B(r)=A- 1 (r) = l- — +1 r-±r 2 , (2) 
r 3 

where M is the central mass, A and 7 are constants. The accepted numerical 
value from current cosmology is A = 1.29 x 10~ 56 cm~ 2 . However, there seems 
to be some ambiguity about the sign and magnitude of 7. Mannheim and 
Kazanas fix it from fiat rotation curve data to be positive and of the order of 
the inverse Hubble length, while Pireaux (2004) argues for I7I ~ 10~ 33 cm -1 . 
Edery & Paranjape (1998) obtained a negative value from the gravitational 
time delay by galactic clusters. We shall keep the value of 7 open, but for 
purely illustrative purposes, take the value 7 = 3.06 x 10~ 30 cm -1 (Mannheim 
2006). We emphasize that the methods adopted here do not need to assume any 
particular value of 7— it is essentially kept free to be fixed by more accurate 
observations. 

In the null geodesic equation, A cancels out giving 



d 2 u „ , , , 7 

—r = -u + 3Mu 2 - 1 
dip 1 2 



u + 3Mu 2 -^. (3) 



We have recently solved the light ray equation using the Rindler-Ishak method 
(Rindler & Ishak 2007), and have shown how 7 and A get mixed up in the deflec- 
tion at higher order (Bhattacharya et al. 2010). Usual perturbative expansion 
up to order M 3 gives the final solution of Eq.(3) as: 

1 sin ip 7 
U =r = — -2 

+ [6 + 3i?V - 3Rj(n - 2<p) cos ^ + 2 cos 2ip - 6.R7 sin ip] 

3M 2 

- [96i?7 + 24R 3 7 3 - 10(2 + 3i?V)(7r - 2tp) cos ip + 32 R-f cos 2tp 

- 20 sin ip - 30R 2 -f 2 sin <p + 3n 2 R 2 -f 2 sin <p - 12TrR 2 -f 2 ip sin ip 
+ 12R 2 j 2 ip 2 sin p - 8ttR-j sin 2ip + 16R-jip sin 2ip + 2 sin 3<p] 

M 3 

+ [9^7(tt - 2ip){R 2 j 2 (n 2 - 150 - 4tt^ + Aip 2 ) - 260} cos</? 

128rt 

+ 2{816 + 2916i? 2 7 2 - 36i? 2 7 2 (7r 2 - 27 - 4irp + 4p 2 ) cos 2<p 

+ 21Rj(tt - 2p) cos3v3 - 4 cos Atp - 111 ORy shop + 540i? 4 7 4 + 356 cos 2ip 

+ 1087r 2 i? 3 7 3 sin p - 360ir Rjtp sin ip - 432irR 3 ~/ 3 ip sin tp + 360i?7</? 2 sin tp 

+ 432R 3 ~f 3 ip 2 sin p - 120tt sin 2ip - 3967ri?V sin 2ip + 90ir 2 R~f sin ip 

+ 792i? 2 7 V sin 2ip - 615R 3 -/ 3 sin ip + 240ip sin 2ip + 126i?7 sin 3^}] , (4) 
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where R is related to the closest distance approach r defined by {tp — tt/2 ) 

1 1 
T*o _ 16i? 3 



[AMR{A - 6i? 7 + 3R A Y) - 8R 2 (Rl - 2) 



+ 



M 2 (33 - 96i?7 + 45i? 2 7 2 - 36i? 3 7 3 ) 
3M 3 



64i? 4 



[152 - 432i?7 + 648i? 2 7 2 - 225i?V + 180i? 4 7 



(5) 



Note that Eq.(4) is a more general solution involving 7 and we can recover the 
SdS vacuole case putting 7 = 0. The Rindlcr-Ishak method requires another 
function A(r, p) = which yields for the present solution 



A(r,<p) = (-r 2 )x 



COS If 

[ 32R 3 



(32R 2 - 3M 2 {20 + 3i? 2 7 2 (10 + (tt - 2ip) 2 } 



+ 32Mi?(9M7 - 2) sin tp) - 6M{3M cos 3tp - 8MRj(ir - 2tp) cos 2tp 

+ (10M - 4i? 2 7 + 9MR 2 -/ 2 )(n - 2tp) sin tp}} 

+ M 3 [18i? 7 {130 + tt 2 (10 + 9i?V) - 4tt^(10 + 9i? 2 7 2 ) + AOtp 2 

+ 3i?V(25 + 12p 2 )} cos^ - 48(10 + 27.R 2 7 2 )U - 2( fi) cos 2 P 

+ 648i?7 cos 3tp + 1620i?77r sin tp + 486irR 3 -/ 3 sin tp - 9ir 3 R 3 ~/ 3 sin ip 

+ 3240i?7<^ sin tp - 972R 3 j 3 p sin tp + 54n 2 R 3 j 3 psm tp 

- l08irR 3 ~f 3 tp 2 sin p + 72R 3 ~f 3 tp 3 sin p - 944 sin 2p - 230AR 2 -/ 2 sin 2p 
+ 1447r 2 i? 2 7 2 sin 2p - 576irR 2 ~/ 2 psm 2tp + 576R 2 j 2 tp 2 sin 2p 

- 1627ri?7 sin 3^ + 324^^ sin 3<^ + 32 sin 4y>]. (6) 

Assume a small angle tp b at the vacuole boundary radius r — r b such that 
s'm tpi, ~ tpi,, and cospb ~ 1. Then the above gives 



n 2 R 



+ M 



2 3tt7 3 7 2 

15tt I27 457T 7 2 97 3 

8^3 ~ "R 2 " + 16i? 4~ 

457 2 93f, 97r 2 7 2 (^ h 



39 



3irj(p b 



16i? 3 2i? 2 

9^7 2 ^ 2 9 7 V 3 ~ 



16R 



+ M 3 
9ir 3 r 



73 

4R 4 



32i? 
11437T7 



+ 



i? 2 

243 7 2 



+ 



+ 



128i? 
2i? 4 
2i? 2 



64i? 3 4i? 2 

^ 135 7 4 157T^, 



8i? 

9^ 2 7 2 

16i? 2 



8i? 

675tt7 3 
64i? 



16 



4i? 4 



7A7jtp b 45ir 2 ~ftp b 



32R 3 



32R 3 



81n^ 2 (p b 6757 3 (^ fc 8l7r 2 7 3 ^ fc 



8i? 2 ' 64i? 
189tt7 3 (^ 2 457^ 



32i? 



8i? 3 



64i? 
997Vg 
16i? 



457T7(p 2 
8i? 3 



(7) 
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and Eq.(6) gives 




-97r 2 i?V - 192i?7^ b + 72TrR 2 -f 2 ip b - 108i? 2 7 2 v? 2 ] + 0(M 3 ), 



(8) 



where the terms 0(M 3 ) are straightforward but rather lengthy and hence not 
displayed here. 

Note that observations give us values of M and R for a lens, but we have 
only one equation (7) for two unknowns ipb and r&. Hence we need to specify 
any one of them from independent considerations. Along with Ishak et al. 
(2008), we shall employ Einstein- Strauss prescription (Einstein & Strauss 1945; 
Schucking 1954) to determine assuming that the vacuole has been matched 
to an expanding FRW universe via the Sen-Lanczos-Darmois-Israel junction 
conditions (Sen 1924; Lanczos 1924; Darmois 1927; Israel 1966). In general, 
the vacuole radius would also change due to cosmic expansion, but we shall 
consider rj, at that particular instant to of cosmic epoch when the light ray 
happens to pass the point of closest approach to the lens. Thus r\, is determined 
by exploiting the Einstein-Strauss prescription [see Ishak et al. (2008)] 



To achieve exact matching with the exterior FRW universe, the energy density 
p within the vacuole should have a contribution from A besides that of ordinary 

matter p m , that is, p = p m + p\ = ^-(tt m + Q A ), where O m = 87rp m /37J 2 , 
= 8ttpa/3Hq are the matter and dark energy densities in dimcnsionless 
form. Current observations suggest that the universe is spatially flat so that 

p = ^critical = -g^ 1 , which in turn imply that il — O m + fl\ = 1- Type la 
supernova observations yield ft m — 0.27 so that £1a = 0.73 (Riess et al. 1998; 
Perlmutter et al. 1999; Carroll 2001; Page et al. 2003; Peebles & Ratra 2003; 
Spergel et al. 2007). For computational purposes, we shall take the density to 



be p in frw = Pcriticai = ^ = 1.1 x 10 29 ( 75 km/ t/M P c ) gm/cm 3 (Weinberg 



1972) inside the vacuole. A slight deviation from this value would not drastically 
alter our conclusions. Normalizing the scale factor to a(i ) = 1 and dropping 
suffixes, the above prescription translates to 



where M is the lens mass often expressed in units of sun's mass M Q = 1.989 x 
10 33 gm = 1.475 x 10 5 cm. 

We should now solve the cubic Eq.(7) in fb to find three roots designated 
by (pi — ifl(p, M, R,j) where i = 1,2,3. For the SdS vacuole (7 = 0), we 
fortunately get only a single root ipi,, which becomes, using the Einstein-Strauss 
prescription for from Eq.(10), 



n in SdS = a(t)r b in FRW? M S dS = ^r r b 



'b in SdS X Pin FRW- 



(9) 




(10) 



<Pb = 



90nM 3 + 96RM 2 - 16(7rp) 1 /3(6M) 2 / 3 J R 3 
117M 3 - 48MR 2 



(11) 
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Our interest is to express the deflection angle ip in terms of M, R and an as 
yet unspecified p. We shall use the above value of ip b later. The Rindlcr-Ishak 
formula for ip at r = r b is 



, , r by /B(r b ) 

tan ^ = — m — • ( 12 ) 

\ A b\ 

From Eq.(7), one sees that r b contains the radius R which is a real root of 
Eq.(5) For large distances, there is little difference between R and tq. So we 
shall identify R with the Einstein radius where the closest approach distance ro 
appears H We are not using the impact paramter here. 

Our general algorithm for calculation proceeds along the following analytical 
steps: 

(1) Put the expression for r b from Eq.(7) and A b from Eq.(8) into Eq.(12) 
for deflection angle. 

(2) Expand the right hand side of Eq.(12) in first power of 7 in order to 
separate out its contribution from the SdS one. Thus we write formally 

tan ?/> total ~ C(tp b , M, R) + jD(ip b , M, R) (13) 

where C, D are known functions to be expanded in powers of M. For small 
V> total , we decompose 

tan V total =a V'' *" 1 = ^ SdS + ^ MKdS - (14) 

(3) Expand both C and D up to the power M 2 to see the individual contri- 
butions of terms. 

III. SdS vacuole: A effect 

This case corresponds to 7 = and we have to be concerned with only 
C = tan-0 sds ~ ^ sdS . The boundary radius r b of the vacuole from Eq.(7) is 
[step (1)], 



n 



16i? 4 



32MR 2 + l6R s cp b + M 2 (30irR - 39R<p b ) + M 3 (292 - 60tt^ + 120^f) 

(15) 



and Eq.(8) gives 



2Mr 2 ip b M 2 r 2 f 39 15nip b 15cpg 



R R 2 R 3 1 16 

M 3 r 2 / 15tt 25<p b \ 



+ 



R 4 { 4 4 J 



(16) 



2 The half angle of the Einstein ring subtended at the observer is defined as 9e = 2eDi s /D s, 
the suffixes o, I, s in D representing angular diameter distances between observer, lens and the 
source, assuming all to be situated on a "line". We have taken R = Re = D oi tanSg ~ D oi 6E 
for small Be- 
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Note that Eqs.(lO), (11) and (15) are consistent. Putting the values of r& and 
Ab in Eq.(12), assuming a small ?/> sds , and expanding in powers of M we get 
[step (2)] 



SdS 



+ M 2 



- 1 



+128.R 3 



Ai? 2 \ 


<fb + M 


6^ ) 


/5ttA 


13A 


Wl 


16^ 


AR 2 \ 


/ M 




\ R 2 




39 




2048i? 5 



: -pb 



(2RA 
2A 



V3^ 



2 2A^ x /2i?A 
_ _ _ , 

240?ri? - 312R(p b 
128i? 3 



1024i? 5 512.R 5 J 



AR 2 

_ fb 
R 



2 2^| 



0(M 3 ). 



(17) 



3M 



To simplify calculations, we shall now expand tp b of Eq.(ll) with p — 
obtaining 

D O Jt/f 

r2 ' 5 (18) 



R 2M 39 15tt\ 



Using it in i/j SdS , and collecting terms of similiar orders in M, we get 



, SriS ARr b M 



R 2 R 2 A 
2+^- — 



M 2 r 
"R 2 " 



15tt 7i? 3 
4^ 



157Tfl 2 



4i? 5?ri? 2 A 

n 



R 3 A 25ARr b 



16 



24r b 



96 



(19) 



As usual, for small angle, taniySft ~ ip b , so that the deflection e sds for nonzero 
tfb is, by definition (Rindlcr & Ishak 2007) 



e sds = tan(^ sds - y b ) ~ V> SdS - ¥>&■ 



(20) 



Assuming that r& >> i?, R » M, and collecting terms of interest, we get the 
total deflection 

Sris ARr b AM 15nM 2 2MAR 5ttM 2 A 2MR , . 

2e =- ^ + + 3 — + ^r- (21) 



Clearly, the above yields the Ishak et al term — A|rt as we n as the well known 
Schwarzschild terms. Interestingly, identifying the constant R ~ ro from Eq.(5), 
the fourth term, viz., t\ c - = — 2M ^ r ° looks numerically like the same as what 
Sereno (2009) calls the local coupling term between M and A, since the term 
does not depend on the vacuole radius r b (see also the Appendix). Nonetheless, 
it does depend on the particular path via R. Remarkably, we also discover a 



new local coupling term in the second order, viz., t\ 



5ttM 2 A 



coming from 



the fifth term in the second square bracket in Eq.(19). This seems to be a more 
genuine local coupling term because it does not involve the parameter R of the 
light trajectory. However, both the terms contribute repulsively to bending. 
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A certain thing is to be noted here. We might start with the first order 
differential equation already containing A through B{r): 



1 



drV 
dip) 



B(r) 



1 



(22) 



where b is the impact parameter defined as £/e. Then one can define b using 
the closest approach distance r = r , where = 0, which yields from the first 
order Eq.(22) the value of b as 



1 



"I 1/2 



B(ro) 



This gives 



r ~ b 1 



f 



M 

T ' 

M 

' T 



PA 
6 

6 



MAb\ 



(23) 

(24) 
(25) 



Using the values of R = ro and 1/R = 1/ro from above into the expression (), 
we find the relevant terms to add to 



2e 



SdS 



bAr b AM 2MA6 



2MA6 



+ terms in M 2 



(26) 



so that the local coupling term 



2MAfc 



cancels out! 



In our opinion, some caution should be exercised in the interpretation that 
such a local coupling really vanishes due to that cancellation. Note that both 
the original Rindler-Ishak (2007) or vacuole (2008) method do not at all use 
the first order path equation in which A already appears. Their whole package 
consists of the second order differential equation in which A does not appear 
(which led people to believe that it does not hence affect light bending) and the 
definition of the invariant angle ip to capture the effect of A, without needing any 
further ingredients. To use the first order path equation (22) already containing 
A at any stage of the present calculation would mean capturing the effect of A 
twice. We argue that both the trajectory equation (3) with R ~ ro and the 
first order equation (22) with b should not be simultaneously used. Integration 
of the first order equation (22) with b and the present vacuole method should 
be treated as mutually exclusive ways, both separately yielding the bending 
with expected local coupling terms. We obtained not only the Sereno-like local 
coupling term — 2M ^ r ° but also a new term — 57rj ^ A including other new terms, 
the most notable one being — AI ^ rb . 

As argued above, the vacuole method has been exclusively tailored to capture 
the effect of a parameter (like A) that has disappeared from the path equation. 
To justify this statement, we might try to assess the influence of the remaining 
parameter 7, which has not disappeared from the second order equation (3). If 
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we still apply the present vacuole method, it can be verified that the known effect 
— ^ on bending does not appear in <0 MKdS (We omit the details). However, 
we can still use the initial Rindler-Ishak non-vacuole method (2007) to obtain 
exactly this effect as already shown in Bhattacharya et al. (2010). Below we 
shall show that even a reverse prescription in the non-vacuole method yields the 
same effect — 

V. The 7 effect by alternative prescription 

To capture the effect of 7 by means of the invariant angle, we assume here 
for simplicity A — and restrict to first order in M. In principle, this case 
corresponds to C = and we have to be concerned with only 

7 ^ = tan^ MKdS ~V MKdS - (27) 

But as stated above, 7D does not yield the effect — Therefore we proceed 
as follows. 

First, a few words about the original Rindler-Ishak (2007) prescription and 
the alternative one we are going to follow here. In the asymptotically non-flat 
metric, the limit r — > 00 makes no sense. Therefore they prescribed that the 
only intrinsically characterized r value replacing r — > 00 is the one at (p = 0. 
The measurable quantities are the various ip angles that the photon orbit makes 
with successive coordinate planes ip = const. While this assumption is perfectly 
valid giving the desired results, we shall implement the Rindler-Ishak method 
following a reverse prescription^], namely, determining ip ^ occurring at a 
point r = 00 on the null orbit and work out how the parameters 7 and M 
appear in the light bending. 

Next, for the unbound orbits associated with lensing, the distance of clos- 
est approach of a light ray to a galaxy will be further from the center of the 
galaxy than the matter orbiting inside it. Hence our goal here is to calculate 
the deflection angle e = ip — tp in the metric B{r) = 1 — — 7- jr under the 
approximation ^ << 1. To this end, we first determine the value of a small 
nonzero <p lying on the null geodesic at r = 00 using the orbit Eq.(4). For small 
(p, we take sin<^ ~ ip, cos(p ~ 1 and neglect terms of ip 2 and higher. Then, at 
r = 00, the orbit Eq.(4) yields 

= -\ + | + ^[8 + 3i?V - 3i? 7 (7T- 2<p) - QR 1V ] 
W 2 

- + 24i?3 7 3 - 10(2 + 3i?V)(^ - 2<p) 

- Utp - 167ri?7^ - 3QR 2 j 2 ip + 3ir 2 R 2 ] j 2 : tp] . (28) 

3 There seems to be no obvious operative way of actually measuring the azimuthal angle 
on the sky unless one is able to perform the experiment twice, once with the lens in the way, 
and once without (just as is done with gravitational bending of light by the sun). However, 
one can define a distance r on the orbit by using a coordinate system in which the lens is 
at r = 0. Then for an observer far enough beyond the lens, the very fact that light reaches 
the observer entails that the point r = 00 does lie on the geodesic. The form of the geodesic 
then determines at what value of tp 7^ the source lies. This is what we are calling reverse 
prescription here. 
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Solving for ip, we obtain 

8MR(8 - 3irRj + 3i? 2 7 2 ) - 16R 3 j + 6M 2 {5tt(2 + 3i? 2 7 2 ) - 4i? 7 (16 + 3i? 2 7 2 )} 



V = 



M 2 (78 - 487ri? 7 + 90i? 2 7 2 + 97r 2 i? 2 7 2 ) - 32i? 2 

(29) 



Calculating for small ip, we find 



ip ~ tan ip 



2du 
' dip 



du\ 1 2M 7 

dip J V r 2 r 3 r ' 

which goes to as r — > oo since ^ 7^ at the value of </? derived in Eq.(29). 
Thus, V = and the one way deflection then is e = — </?, which easily expands 
to 



2M 



, 15ttM 
1 + 



16R 



7 



i? 3ttM 423M 2 



2 4 32i? 



(31) 



We have checked that this result exactly coincides with that obtained by 
the perturbativc Bodcnner-Will perturbative method (2003). We find that all 
terms in the second square bracket are positive, meaning that the effect of 
7 > is to diminish (and 7 < is to enhance) the Schwarzschild bending even 
up to second order in M. We also find that Eq.(31) nicely reproduces the one 
way deflection = — 2E^j obtained by Edcry and Paranjapc (1998) using 

Weinberg's method. As mentioned, the same result (31) follows also from the 
unaltered Rindler-Ishak prescription (2007) as well. We shall now discuss some 
physical implications of Eq.(31). 
VI. Physical implications 

First note that in the halo we have obtained a new coupling term 
between M and 7, which is independent of the trajectory parameter R. Next, 
the term — shows repulsion for 7 > 0, which is consistent with time delay 
investigations (see e.g., Edery & Paranjape 1998) and attraction if we choose 
7 < 0. We emphasize that we are not concluding anything about the correct 
sign of 7, which must be decided by independent observations. When M = 
and 7 > 0, we obtain a negative (repulsive) bending of light or e = — which 
coincides with the conclusion by Walker (1994). 

On the other hand, in the galactic halo region, where R > i?E (the Einstein 
radius) and ^ « 1, one would like to obtain a positive (attractive) light 
bending there. This is possible only if one assumes the condition 

2M 7 i? 

£ > ° ~R~ > 2 (32) 

to hold. Accurately observed lensing data by galactic clusters are now available. 
We then find from Table I that the observed values of M and Re do indeed 
respect the inequality (32). Clearly, even if pure 7 > leads to repulsion, in 
the competition between this repulsion and Schwarzschild attraction, the latter 
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might win leading to the impression of an overall attractive bending. This can 
happen in the lensing by galactic clusters, as described in the table below. 

Lens data for M, Re and references are taken from Ishak ct al. (2008), and 
converted here to length units using M© = 1.475 x 10 5 cm, 1 kpc = 3.0856 x 
10 21 cm. The references are as follows. 1: Abell 2744 (Smail et al. 1991, Allen 
1998), 2: Abell 1689 (Allen 1998, Limousin 2007), 3: SDSS J1004+4112 (Sharon 
2006),4: 3C 295 (Wold et al. 2002), 5: Abell 2219L (Smail et al. 1995a; 
Allen 1998), 6: AC 114 (Smail ct al. 1995b; Allen 1998). We shall take the 
rotation curve fit value 7 = 3.06 x 10~ 30 cm -1 purely for illustrative purposes 
and A = 1.29 x 10~ 56 cm~ 2 in both the tables below: 

Table I 



Cluster 


M (cm) 


-Re (cm) 


2M 
Re 






7 -Re 
2 






Abell 2744 


2.91 x 10 18 


2.97 x 10 23 


1.95 x 


10 


-5 


0.45 


x 10 


-6 


Abell 1689 


1.38 x 10 18 


4.26 x 10 23 


0.65 x 


10 


-5 


0.65 


x 10 


-6 


SDSS J1004+4112 


6.28 x 10 18 


3.39 x 10 23 


3.70 x 


10 


-5 


0.51 


x 10 


-6 


3C 295 


10.5 x 10 18 


3.94 x 10 23 


5.32 x 


10 


-5 


0.60 


x 10 


-6 


Abell 2219L 


4.75 x 10 18 


2.66 x 10 23 


3.57 x 


10 


-5 


0.40 


x 10 


-6 


AC 114 


1.36 x 10 18 


1.68 x 10 23 


1.61 x 


10- 


-5 


0.25 


x 10 


-6 



It is evident from the above table that the term is smaller than the 
Schwarzschild term so that the overall bending is always attractive for 
7 > 0. One might want to have an idea of the radius R where the leading 
order Schwarzschild and 7— bendings balance each other. The value of i?b may 
be taken to demarcate the boundary of the halo dark matter surrounding each 
individual cluster. This happens at 

[m 

R h = 2J— cm. (33) 

The deflection e below R < R\> is always attractive, as should be the case. 
Table II shows that the halo boundary i?b can be several times larger than 
Re- However, the values of i?b tabulated here rely crucially on the value of 
7 and if its value is lowered by one order of magintude than considered here, 
Rh will increase by that order. Conversely, if one particular halo boundary is 
observationally determined, then it would provide us with a determination of 
7. One could then examine if that new value of 7 explains i?b of other clusters. 
If it does, then it would support Weyl theory. Observations seem as yet far too 
inconclusive about the sizes of the halo. 

Table II 



Cluster 


M (. 


:m) 


-Re (cm) 


Rh ( 


cm) 


Abell 2744 


2.91 


X 


10 18 


2.97 x 


10 23 


19.5 


x 


10 23 


Abell 1689 


1.38 


x 


10 18 


4.26 x 


10 23 


13.4 


x 


1Q 23 


SDSS J1004+4112 


6.28 


x 


10 18 


3.39 x 


1Q 23 


28.6 


x 


1Q 23 


3C 295 


10.5 


X 


10 18 


3.94 x 


10 23 


37.0 


x 


10 23 


Abell 2219L 


4.75 


X 


10 18 


2.66 x 


10 23 


24.9 


x 


10 23 


AC 114 


1.36 


X 


10 18 


1.68 x 


10 23 


13.3 


X 


10 23 
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Though our interest so far has only been in the galactic clusters, one might 
still want to compare the magnitudes of the 7— related effects with the Schwarzschild 
ones in the solar system although the region around the Sun contains galactic 
matter. For a light ray grazing the Sun, we have the following numerical values: 

M Q = 1.475 x 10 5 cm, R Q = 6.96 x 10 10 cm, 7 = 3.06 x KT^cnT 1 (34) 

so that 

— ° - 4.24 x 10- 6 , — = 2.65 x KT 11 , (35) 

tlQ ID-Kq 

2^ = 1.06 x 10" 19 , = 1.06 x 10 - 24 . (36) 

2 4 v ' 

We find that the 7— correction terms are considerably small compared to ^rj 2 -, 
therefore the effect of 7 would be negligible near the Sun. However, as illustrated 
in Table I, the effect of 7 near any galactic cluster scale is not as negligible. The 
fact that 7 is meaningful only on such large scales has been conjectured in 
the literature, but here we have found its support from a completely different 
viewpoint, viz., from the Rindlcr-Ishak bending. 
VII. Summary and results 

We calculated light deflection in the vacuole model up to third order in 
M and confirmed that the extension of Rindler-Ishak method to the Einstein- 
Strauss vacuole, as originally developed by Ishak et al. (2008), reproduces the 
Schwarzschild M— dependent bending terms as well as the A — dependent 
terms, see Eq.(21). In particular, we have found a local coupling term — 2M ^ r ° 
similar to that by Sereno. We have also found a more interesting coupling term 
_5ttM_a mc i uc i m g other new terms, the most notable one being — A ^ rb . It 
would be of interest to discuss the recessional impact too (Ishak & Rindler, 
2010), but it requires a separate and detailed investigation. 

The idea of a cut-off transition region between the halo boundary and the ex- 
terior dS cosmology was conjectured, but not implemented, by Edery & Paran- 
jape (1998) over a decade ago. The SdS vacuole model by Ishak et al. (2008) is 
philosophically the same in idea but different in content. It envisages a transition 
radius r& between the SdS vacuole boundary and the exterior FRW cosmology 
implementing the Einstein- Strauss suggestion. The vacuole surrounding the 
lens should be deviod of matter, and therefore the model particularly applies to 
galactic clusters rathen than local objects like the Sun, which is surrounded by 
galactic matter. 

We have argued that the vacuole method is exclusive to cases where the cos- 
mological constant A disappears from the second order differential path equa- 
tion. To exemplify it, we applied the vacuole model in the calculation of the 7— 
dependent effects in Weyl gravity. We note that the parameter 7 does not dis- 
appear from the path equation, and thus the vacuole method does not yield the 
known Weyl term — To this end, we point out that the earlier Rindler-Ishak 
(2007) prescription in their non-vacuole method did nicely yield the otherwise 
known Weyl term (Bhattacharya et al. 2010). In the present paper, we showed 
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that an alternative prescription on the azimuthal angle lying on the null orbit 
also reveal the influence of the Schwarzschild (M) and conformal sector (7) on 
light deflection [See Eq. (31)]: It reproduced the correct Schwarzschild bending 
terms due to M > as well as those due to the conformal Weyl parameter 
7. In particular, the known term — 7^ followed exactly. Also we have found 
a new local coupling term 37r ^ between M and 7, which is independent of 
the trajectory parameter R. We chose (not mandatorily) the value obtained 
by Mannheim (2006) from the fit of the galactic flat rotation curve data and 
applied it to the accurately observed data on several galactic clusters taken from 
Ishak et al. (2008). We have shown in Table I that, for R E < R < Rb, the 
light bending is attractive since e ^= ^jf — is always positive masking the 

purely negative Weyl 7— term, while Table II gives possible sizes i?b of the halo 
if the chosen value of 7 is relied upon. Although galactic halo can be modelled 
in many ways [see, for instance, the brane world model, Nandi et al (2009)], 
the interpretations of Weyl gravity in this regime seem as yet conclusive, to our 
knowledge. 
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Appendix 

The integration of the first order equation reads 

^ SCrCn ° - ± / % [I, + | - ^ + ™] • - ± / W> A, b, r)dr (Al) 

It can't be integrated in a closed form. So expanding the integrand / in first 
power of M, we have 

AV 1/2 M ( 1 1 A x ~ 3/2 



f = ^{ ¥ -^~z) -prU-^-sJ = h + h (sayMA2) 

Then, to first power of A, 
h = I fidr 

y/b 2 (3 + r 2 A) -3r In r - In 2 - In {by/3 + ^(3 + r 2 A) - 3r}] 





V3r 




b 3 


3b 5 Ab 3 


Ab 5 


6r 3 


40r 5 6r 


12r 3 


h 


= [ f2dr 





b 2 Ab 2 
r 7 3 

^ + imaginary terms. (A3) 



M[6 2 (3 + 2r 2 A) - 6r 2 ] 



3br*yjl - £ - 



Ab 2 
3 



2M MAb Mb 3 Mb 5 MAb 5 

~b 3~~ + + ~W + ~8^' ^ ' 
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Collecting real terms, we get 

Scrcno _ 2M b MAb b 3 Mb 3 36 5 Ab 3 



b r 3 6r 3 4r 4 40r 5 6r 
MA 2 b 3 Ab 5 Mb 5 MAb 5 



+ -^r + -sis-. (A5) 



36 12r 3 4r 6 8r 4 
which seem to yield that the local coupling term is — 2MAb 



15 



